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Abstract The electron spectrum in the normal phase of cuprates is analyzed 
in terms of the boson-fermion model. It is argued that the existence of the un- 
condensed pairs and the quasi-two-dimensional nature of the crystal structure 
are responsible for the pseudo-gap above Tc observed recently. 

1. Introduction 

Recently, ARPES experiments revealed a novel property of the 
underdoped cuprate superconductors in their normal phase. There 
is a pseudo-gap in the electron spectrum, i.e. a suppression of the 
spectral weight near the Fermi level above the transition temperature 
[1] . It is speculated by a number of authors that the phenomena may 
be attributed to the existence of the uncondensed pairs in the normal 
phase. The present article is an attempt in this direction. 

A common feature of the cuprate superconductors is that the co- 
herence length is comparable with the lattice spacing on the Cu02 
plane. Therefore the Cooper pairs, if they exist, are highly localized 
in the coordinate space and consequently can be regarded as boson 
degrees of freedom. The superconductivity is thereby associated with 
the kinematical Bose-Einstein condensation of these pairing states and 
there must be uncondensed pairs in the normal phase. 

The crystal structure of all cuprate materials is highly anisotropic 
It consists of parallel layers of Cu02, the minimum distance between 
oxygens (where the doped charge carriers may reside) on the same 
layer is less than 3A. The separations between neighboring Cu02 lay- 
ers are not even but periodic with an average spacing around QA. The 
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number of layers in each period varies from one for La2-xSrxCu04 
{Tc < 40K) to three for Tl2Ca2Ba2Cu2,Ox {Tc < 120K). Although 
the coherence perpendicular to the Cu02 layer is required for the con- 
densate to form below Tc, the normal phase properties of the under- 
doped materials are largely two dimensional in nature. This has been 
confirmed by a number of experiments on the transport coefficients 
of cuprates [2]. 

The Bose-Einstein condensation is characterized by the relation 
that the thermal wavelength of the bosons, 



At^J^ (1.1) 

with the effective mass of bosons, at the transition temperature 
T — Tc should be comparable with the inter-boson distance, I, so 
that quantum mechanical coherence takes place. The ratio ^ is 1.38 
for an ideal Bose gas and 1.65 for Hell at the A-transition. For high 
Tc materials, the ratio 



I 



2.8 (1.2) 



was estimated in [3] based on the results of the fiSR experiment [4]. 
The reasoning went as follows: According to the experimental results, 
the transition temperature Tc of the under-doped material is inversely 
proportional to the square of the magnetic penetration depth with a 
constant of proportionality universal for all cuprates. If the supercur- 
rent is identified with the superfiuid of bosons. A relation 

Tc{K) = 10^!^^ (1.3) 

was obtained, where rUe is the electron mass in vacuum, m^, is the 
in- plane effective mass of the bosons, nt is the density of the bosons 
(counted as number of bosons per and per layer) and d the average 
separation between Cu02 layers (in the unit of A). Combining the 
relation (1.3) with d = 6 A, ni, = and the definition of the thermal 
wavelength (1.1), the ratio (1.2) follows independent of nif,. This ratio 
is certainly consistent with the picture of Bose-Einstein condensation. 

2 

In an isotropic situation, a relation Ay,-. ~ I implies Tc oc in 
accordance with (1.1) and ni, = l^^. This is not the case of (1.3). 
Therefore we have to resort to the quasi-two dimensional layer struc- 
ture of the system. Consider an ideal Bose gas in quasi-two dimensions 
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whose kinetic energy consists a free motion term in x — y plane and a 
hopping (tunneling) term in 2;-direction, i.e. 

^p,K = ^ + 2^6(1 - COS Kd), (1.4) 

where p is the momentum in x — y plane, K is the Bloch momentum in 
2;-direction and ti, is the out-of-plane hopping amplitude. Expanding 
the cosine of (1.4) to term, we obtain the out-of-plane effective 
mass 

Mh = (2^2)-! (1.5) 
The density of states of the spectrum (1.4) is 



rub 



1 for a; > 4tft; 



P(^) = 1 2 .^-1 r . (1.6) 



271 \|sin '^^^ foTU<4tb 

The boson density (counted as number of bosons per unit area and 
per layer) at a temperature T and a chemical potential < is 

'''~J-.2^J (27r)2e/3(-M)-l - A2:y_,2^ 1 



g-/3[2tb(l-cos6>)-/u]' 

(1.7) 

where P = and Ay = -y/ ^^^j, is the thermal wavelength of the 

bosons in the x — y plane. For f3ti, << 1 and -~(3fi << 1, the integral 
(1.7) may be expressed in terms of elementary functions. We find that 

1 , 2kT , , 

n?, = In , - (1.8) 

(2tfe-yu) + ^(2t,,-/i)2-4t2 

At the transition = and we have 

Regarding x — y plane as a Cu02 layer, it was found in [5] that 
this quasi-linear relation between and Tc fits the //S'i? result (1.3) 
remarkably well with 

In ^ ~ 8 at Tc = lOOK , (1.10) 



which imphes that 



tb/nTc = 3.6 X 10 



-4 



(1.11) 



Furthermore, the formula (1.10) can be generalized to the material 
with multi-layers in a period in z direction provide we replace the hop- 
ping amplitude ti, in (1.10) by the geometrical mean of the hopping 
amplitudes between each pair of neighboring layers within one period. 
This also explains the universality of the n^, — Tc profile (1.3), since 
the tunneling nature of boson motion in ^-direction makes the geo- 
metrical mean of the hopping amplitudes sensitive only to the average 
separation between the Cu02 layers and this average is approximately 
the same for different cuprates. 

The large logarithmic factor is obviously caused by the large fluc- 
tuations of low-lying bosons, i.e. the 2D nature of the density of 
state (1.6) down to the tiny energy scale tb. Since the square root 
of the logarithm of (1.10) gives the quoted ratio (1.2), the quasi-two- 
dimensional character also responsible for the factor two enlargement 
of this ratio compared with that of an ideal Bose gas in three dimen- 
sions, bearing in mind that within 3D, the change of this ratio from 
an ideal Bose gas to the highly interaction Hell is merely 20%. 

A simple phenomenological model dealing with system with local- 
ized Cooper pairs was proposed in [3] and [6], which will be refered 
to as the boson-fermion model. In the rest of the paper, we shall first 
review this model and then apply it to explore the consequences of 
the presence of uncondcnsed bosons within a quasi-two-dimensional 
crystal in the normal phase. 

2. The Boson-Fermion Model in Quasi-Two Dimensions 

The boson-fermion model appropriate to the quasi-two dimen- 
sional crystal structure is formulated in ref.[7]. The grand Hamilto- 
nian of the system reads 



where ap^ and ap {bp and bp) are annihilation and creation operators 
of electrons (bosons), the subscript s denotes the spin orientation 




P,s 



P 
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(s =T or i), 

€0= :^ + 2tf(l- cos Kd)- u, (2.2) 
^ zruf ■' 

2 

uj0^ ^ + 2tb(l-cosKd) + 2(u- ij), (2.3) 
^ 2mb 

mf{mi)) is the in-plane effective mass of electrons (bosons), tf{ti,) is the 
out-of- plane hopping amplitude of fermions (bosons), is the chemical 
potential of the system, 2i/ is the energy of a static boson relative to 
two static electrons {2u = (jJp\ p^^ — '2ep\ p^g), ^ is the total area of a 
Cu02 layer and M is the total number of layers of the whole system. 
The vectors with capital symbols, P and Q are all three dimensional, 
i.e. 

P = p + Kz 

and 

Q = q + K'z 

with q parallel to the x — y plane, K' Bloch momenta in z 
direction whose value varies between —iz/d and ir/d. The conserved 
electric charge is given by 

P,s P 

For resonant bosons, > and the dimensionless coupling constant 
g is related to the boson width in vacuum through 

g2 (I iov V > 4:tf; 

T = -^ruf I I g-j^-i otherwise. 

For 5^ = 0, the model implies an ideal Bose gas and an ideal fermion 
gas in chemical equilibrium. In what follows we shall assume that 
V > 4tf. The effective expansion parameter of the perturbation series 
is g defined by 

f-l-.^^ (2.5) 

The condensation mechanism with resonance bosons can be un- 
derstood easily in the weak coupling limit, g « 1 and at T = 0. 
Starting with low charge density, the fermionic levels ep are filled 
first. The bosons are unstable, existing in virtual state through the 
interaction and the condensate is of BCS type. With the increasing 
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charge density, the fermion levels are filled up to y. Additional charges 
will be stabilized at the zero momentum boson level to minimize the 
energy since the decay channels of a static bosons are all locked by 
the Fermi sea and the condensate is of Bose-Einstein type. 

In the super-phase, we introduce a long range order B representing 
the Bose condensate at the P = 0, i.e. 

feplp^o = ^^e*" + 6', (2.6) 
where S > 0, a is a constant phase, and then 

H = Ho + Hi, (2.7) 

P 

+ X^[ep(4T«FT + + A(e'''4T''-Pi + e"'''«-n«pp] (2.8) 

p 

with A = gB, and Hi includes the cubic terms in a's and Vs. The 
electron spectrum to the zeroth order in Hi but to all orders in A 
can be obtained via a Bogoliubov transformation of Hq. In order 
to compare notes with the normal phase results later on, we adapt 
diagrammatic technique to derive the electron spectrum. All the di- 
agrams throughout this paper are thermal with discrete imaginary 
energies (Matsubara variables) and continuous spatial momenta. The 
real-time (real energy) Green's functions can be obtained by replacing 
all external Matsubara energies by real and continuous energies. The 
uniqueness of such a replacement was established in [8]. 

Returning to i^o, we associate a solid line with the bare electron 
propagator 

(2.9) 

IVn - ep 

with = '^{'^ \) and ep given by (2.2), a cross with two incom- 
ing fermion lines is associated with a factor — zA and that with two 
outgoing lines with a factor zA. The electron propagator to all or- 
ders in A but to zeroth order in Hi (represented by a solid line with 
two arrows), with a momentum P and a discrete energy Sn{P), 
corresponds to the sum of the diagrams in Fig.l. 
We have 

S„(P) = -J— f (_iA)"(iA)"(^^)"(^)" 
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+ » X « X » X « x » + ••• 

Figure 1. The diagrammatic expansion of the electron propagator in the super 
phase 

= 2 J (2-10) 



where Ep = + A^. Replacing the Matsubara energy iun hy 
Pq + iO'^ with pQ real, we obtain the retarded electron propagator 

The corresponding spectral density is 
A{po,P) = -ReSR{po,P) = ^-^^-^[d{po-Ep)+S{po+Ep)] (2.12) 
with two infinitely sharp quasi-particle peaks. On writing 

-5(^0, P) = ' ^ (2.13) 

po + iO+ -€p-i:{pQ + iO+,P) 

we extract the self-energy function 

S(po,P) = ^^ (2.14) 



3. The Electron Self-Energy in the Normal Phase 



Above the transition temperature, T > Tq, the coherence between 
different layers is no longer crucial and we shall approximate the model 
Hamiltonian (2.1) by a purely two-dimensional one, i.e. 

p,s p 
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where 



9_ 

p,q 

2 



6,^=^-^, (3.2) 

u^P = £^^+S (3.3) 

with 6 — 2{iy — n) whose value in the vicinity above Tc is regarded 
of the same order of th in (1.11). All momenta are two-dimensional 
from now on. 

3.1 The self-energy to one-loop order 

In terms of the thermal diagrams, a free electron propagator, rep- 
resented by a solid line, with a momentum p and an energy z^(n+ |) 
corresponds to 

with n an integers; a free boson propagator, represented by a dashed 
line, with a momentum p and an energy i^n corresponds to 

— ; (3.5) 

t-^n — ujjj 

a vertex with two incoming electrons and one outgoing boson corre- 
sponds to —ig and that with two outgoing electrons and one incoming 
boson to ig. The energy and momentum are conserved at each vertex. 
A loop corresponds to a sum and an integral 



P^J (27r)2 

with I the integer labeling the Matsubara energy of a chosen line of 
the loop and q the momentum of that line. A minus sign is associated 
with a fermionic loop. All arrows in the diagrams follow the charge 
flow. 

The electron self-energy to the one-loop order is given by i times 
the amputated part of the diagram in Fig. 2a, i.e. 



^n{p) = i^{-ig)igY^ 



r d^q 


i 


i 


1 (27r)2 




U^(Z-n-i)-w^-_pJ 
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(a) (b) 

Figure 2. The one-loop diagrams for the self-energy of an electron (a) and the 
self-energy of a boson (b) . 

where N\) and Nf are boson and fermion distribution functions given 

by 

and 

NfiP) = (3-8) 

e'^ p + 1 

respectively. With the Matsubara energy i^(n + |) replaced by po + 
iO'^, we obtain the retarded self-energy function 

^ 2 f d\ Nt{q)+Nf{q-p) 

^{Po,p)=9 77r-yl — ^n^- (3-9) 

J (27r)^po -c^g - e9-p + *0+ 

If the bosonic chemical potential —5 were zero, the integration over q 
would be logarithmically divergent. Therefore, in the limit 5 — > 0, we 
expect 

S(po,p) ^ (3.10) 

with 

A2 = ^!^ln^. (3.11) 

The right hand side of (3.10) is of the same form as (2.14) and a 
perfect gap emerges. With a moderate logarithmic factor, the gap 
becomes a pseudo one as smeared by the imaginary part of S(po,p)- 
The integration in (3.9) can be carried out analytically in the situation 
ji» kT and ji» and ji » po. The result reads 

2 r kT 



^4:ri^6+{pQ + e^-6y 



sign(po + ep - 5) 



X In 



— ITT 



(3.12) 



(3.13) 



r = mi,/mf and ^(1/2) = —1.4604. For p on the Fermi surface, 
p = pfi ep'= and we find tfiat S(pO)P)|p=pF ~ '^(Po) + ^'^^(Po) where 



where the function / is defined as 

9 roo 

m = 



ze 



IT Jq [1 + ze ' 



u{po) = 9 



2^ 

27r 



reTsign(po - ^) A/4r/i5 + (po - (5)^ + \PQ - ^\ 



^Arjid + (po - f^)^ A/4ryU(5 + (po - f^)^ - \po - 5| 



i^+l;(e-/^(Po-.)). 



rfJL 



(3.14) 



and 



with 



^(po) = -g 



2^ 

27r 



tikT 



^J^r^^ + (po - 5)2 



+ l(/(e/^(Po-.)) + C(^ 



The spectral function at the Fermi surface is 

1 



TT PO - w(Po) - ^^'(Po) + w((^) ' 



(3.15) 



(3.16) 



where we have renormahzed the chemical potential in (3.16) by sub- 
tracting the constant term u{S). The function ^(po,p) on the Fermi 
surface is plotted in Fig. 3 for In^ = 8, 6 and 4 respectively. The 
opening of a pseudo gap with increasing lnAcT/5 is clearly shown. 

3.2 A resummation of higher order diagrams 

The validity of the perturbative results for E(po,p) requires both 
g^ << 1 and g'^ln^ << 1. Nonperturbative effects will enter the 
game if the second inequality fails. In what follows, we shall examine 
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Figure 3 The spectral function at = 0, with = 10, r = 2 and g'^m},j3/2-K = 1. 
The sohd line corresponds to In kT/S = 8, the dashed line to In kT/S = 6 and the 
dotted line to ln«;T/^ = 4, where x = Po/kT and y = A{po,p)\p^ppKT. 

the higher order diagrams and derive an expression which can be 
extended to the region g'^ << 1 but In ^ ~ 1 

A diagrams of an electron propagator to an arbitrary order con- 
sists of at least a string of electron lines, refered to as the main string, 
joining the initial and the final states. A number of vertices has to be 
attached to this line. In what follows, a vertex with an outgoing boson 
line will be refered to as a source and that with an incoming boson 
line will be refered to as a sink. On account of the charge conser- 
vation, there should be an even number of vertices with sources and 
sinks alternatively attached to the main string. Without decorating 
the boson lines, there are no electron lines outside the main string. 
Unlike the super phase (Fig. 1), in which the boson lines all termi- 
nated at the condensate, a boson line emitted from a source can land 
at any sink of the diagram. Therefore there are N\ such diagrams to 
the 2A^-th order in g and loops in each diagram. Such an expansion 
to 6th order in g is displayed in Fig. 4. Each loop contains a boson 
leaving the main string and then coming back. The Matsubara ener- 
gies of these bosons are denoted by i^lm (m = 1, 2, AT) with Im in- 
tegers and the corresponding momenta by Qm- Im^s are to be summed 
and Qm^s are to be integrated. The leading divergence as (5 — comes 
from the integrals over qm's of the term with /i = /2 = ... = /at = 0, 
where each momentum integral over qm is logarithmically divergent. 
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+ • • • 

Figure 4. The diagrammatic expansion of the electron propagator in the normal 
phase. 

For a tiny 5, each integral remains dominated by the lower limit, and 
contribute the amplitude of the diagram a factor In y >> 1 with A 
a cutoff energy ~ kT. To the leading order of these logarithms, the 
terms with some of l^s nonzero can be dropped and the deviation 
of the momentum of an internal electron line along the main string 
from the external momentum p due to the boson emissions and ab- 
sorptions can be neglected. Collecting all the fermionic propagators 
and vertices with such an approximation, we find that all 2A^-th order 
diagrams contribute equally to the full electron propagator. Snip)-, a 
term 




(3.17) 
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where i^n = '^{n + \) and 



X2 2^bl^T A 2^bVKT A 

A = ^ — In - = 25r In - 

ZTT rUf TT 



with g defined by (2.5). On summing over N , we obtain 
where function F is given by the asymptotic series 

m - (3-19) 

N 

which is formally Borel summable [9] (It is advantageous in this case 
not to distinguish the one particle irreducible diagrams with the re- 
ducible ones). On substituting 

N\= dtt^e-^ (3.20) 
Jo 

into (3.19) and interchanging the order of the integration and the 
summation, we find 

r°° e"* 11/ 1\ 

^«) = /„ *rTa = -r'K-c)- (^'^i' 

where the function Ei is the exponential integral defined in [10] 

Ei(^) = / dt— (3.22) 

J —oo t 



with arg(— 2;) < tt. Ei{z) is analytic on the 2;-plane cut along the 
positive real axis with the discontinuity across the cut given by 

Ei(a; + iO+) - Ei{x - iO+) = -27ri. (3.23) 

The asymptotic expansion with a large \z\ or a small \z\ are given 
respectively by 

T Sn=0 ^ for \z\ » 1 
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with 7 the Euler constant. 

Replacing the imaginary energy ii/^ by pq + iO'^, we obtain the 
retarded electron propagator: 



Sr{po,P)=i ^2 a2 Ei(^ ^2 )^ (3-25) 

which returns to the free propagator in the limit A — > 0, i.e. 

lim Sr{po,p) = \ .^,1 • (3.26) 

The self-energy function extracted from (3.22) according to (2.13) 
reads 

^{P0,P) = \z - z-^ lnEi(^)l (3.27) 



dz 



with z = -^2^, which agrees with the one- loop result if << 
^^In^ << 1. 

For a fixed p, Sji{pQ,p) is an analytic function on the po plane with 
two cuts running from le^-l to oo and from — \ejj\ to — oo. The spectral 
function is given by the discontinuity across the cuts. It follows from 
(3.23) and (3.25) that 

A{Po,P} = -Re5i?(po + iO^,p) 

0, for \po\ < \e^\; 

^5^^e a2 , otherwise. 

The dependence of A{po,p) of (3.28) on the coupling constant g in 
nonperturbative. On the Fermi surface, = 0, we have 



\Po\ 



„2 



AipQ,p)\p=p^ = ^e ^. (3.29) 

Two peaks at po = ±A/-\/2 together with a depletion of states at the 
Fermi level po = corresponding to a pseudo-gap emerge. 
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4. Comments 

In this final section, we shall comments on the validity of our 
results for a realistic system and its relation with known results in the 
literature. 

The electron self-energy function in the boson-fermion model has 
been investigated in [11] and a pseudo-gap was found in D = 1 and 
2 by numerical solutions of the self-consistent equations for T,{po,p). 
The approach taken in this article is purely analytical with emphasis 
on the quasi-2D nature of the crystal structure. The one-loop result 
shows the trend of a pseudo-gap when the logarithmic factor In j be- 
comes sizable and the resummation of the diagrams in Fig. 4 gives 
the nonperturbative effect when the logarithmic factor becomes dom- 
inant. We notice that the entangled diagrams, e.g. the last one of Fig. 
4, which is often neglected in the self-consistent method, contribute 
equally as the others to the leading order of the logarithm. Neverthe- 
less, the result as to the existence of a pseudo-gap of the present work 
agrees with that of [11] qualitatively. 

The boson-fermion model is a phenomenological one, in which 
bosons are regarded as elementary. In a realistic system, bosons cor- 
respond to a pair of electrons and do not represent new degrees of 
freedom. This is not not a serious problem as long as the charge 
density is not too high to cause substantial overlapping among the 
pairs. A theorem proved in [12] states that for any fermionic Hamil- 
tonian, there exist an equivalent boson-fermion Hamiltonian with the 
elementary bosons of the latter correspond to the fermion pairs of 
the former. For a Hubbard like model with an on-site attraction and 
a nearest neighbor repulsion, the low density behavior is indeed de- 
scribed by the boson-fermion model. 

Although the motion of electrons and pairs in the x — y plane 
is assumed to be in continuum, the result will not be changed if we 
replace the motion m. x — y plane by lattice hopping. Furthermore, 
the model Hamiltonian (2.1) can also accommodate a momentum de- 
pendent coupling say 

gocpl-pl, (4.1) 

which results a pseudo-gap with d-wave character. 

In the perturbative expansion (3.17), all boson propagators are 
kept undecorated. The decorations come in two ways: One is to 
replace the bare boson propagators by the dressed ones and the other 
is to include the scattering diagrams of these bosons. Both require 
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fermion loops in addition to the main string of fermion lines under 
discussion. 

The dressing of the boson propagator is through the self-energy 
function Unip) defined at momentum p and Matsubara energy i^^. 
The one loop diagram for n„(p) is depicted in Fig. 2b. We find, to 
that order and n = that 



as p — > with a and b functions of T and /x. The result is believed to be 
true in multi- loops. Since the contribution to (3.14) comes only from 
the boson propagator with zero Matsubara energy and low momen- 
tum, the efi'ect of the boson self- energy is merely a renormalization. 
The results of the last section is not modified if we regard and 5 
being renormalized quantities. 

The diagrams with boson scattering look problematic and the one 
of such diagrams to the lowest order in g is shown in Fig. 5a. 
A naive power counting indicates a linear divergence in the limit S 
0. This however can be eliminated on summing over the diagrams with 
repeated scatterings between two bosons as is shown in Fig. 5b-c. Our 
result will survive if this is also true for multi-boson scatterings. 

The Coulomb interaction has been left out throughout the dis- 
cussions. Strictly speaking, the Coulomb interaction is perturbative 
only if the inter-particle distance is short compared with the Bohr 
radiu [13]. Therefore, the results of the preceding section would be 
unaffected in a hypothetical situation when the charge density is high 
enough and the bosons are truly elementary. The situation with the 
real materials is, however, not in the perturbative region at all. On 
the other hand, the shielding of the Fermi sea truncates the long range 
tail of the direct Coulomb force and there are possibilities that our 
results survive. Diagrammatically, inclusion of the Coulomb interac- 
tion amounts to decorate the electron and fermion propagators and 
the boson-fermion vertices and to add scattering diagrams between 
electron and bosons. As long as the Fermi-liquid behavior of electron 
spectrum and the low momentum behavior of the boson spectrum 
are not modified by Coulomb interactions alone, our result remain 
qualitatively intact. 

Technically, the most drastic step leading to (3.25) is the Borel 
summation of the perturbative series (3.18). Is the sum unique? An 
asymptotic expansion of a given function -^(C), i-e. 



Un{p) = a + bp'^ 



(4.2) 




(4.3) 



n 
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(c) 

Figure 5. (a) The lowest order diagram of the electron propagator with scattering 

of virtual bosons; (b) The propagator diagram with the scattering between the 
two virtual bosons to all orders; (c) The expansion of the shaded part of (b). A 
cross-shaded bubble is a two particle irreducible part. 

may well be the asymptotic expansion of another function, say -F(C) + 
_ 1 

e ^ . Therefore the sum of an asymptotic series is in general nonunique 
unless certain conditions are preimposed on the sum. According to 
the Watson theorem in [9] , the sufficient conditions for the Borel sum 
to be a unique function F{() which produce the asymptotic expansion 
of the right hand side of (4.3) are: 1) F{Q has to be analytic within 
the region D defined by \(\ < r and |arg(0| < ^ + a with r > and 
a > 0; 2) \an\ = 0{n\a''y, with a > 0; 3) 

N 

\F{0 - E = 0{{n + l)!(7-+V-+i). (4.4) 

n=0 
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uniformly within D. It follows from the rigorous spectral representa- 
tion of the fermion propagator that the condition 1) is indeed satisfied 
if we identify z with A^/ {i^n + ep) and the coefficients of the pertur- 
bation series (3.18) meet 2). The condition 3) is difficult to justify 
nonperturbatively. Furthermore, each term in (3.18) is only a part 
of the perturbation in terms of g, since only the leading logarithmic 
term is reserved. Therefore the rigorous mathematical justification of 
the Borel sum (3.18) remains open. 
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